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ABSTRACT 


A combination  of  computer  analysis  and  scale  model  testing  was  utilized  to 
develop  a nozzle  which  would  increase  the  performance  of  thrust  augmenting 
ejectors.  Scale  model  tests  were  conducted  on  various  multi-lobed  and 
vortex  generating  nozzles.  Predicted  jet  characteristics  were  obtained 
by  calculating  a finite  difference  solution  of  Reynolds  equations  for  the 
three  dimensional  flow  field.  A two-equation  turbulence  kinetic  energy 
model  was  used  for  closure.  It  is  demonstrated  that  the  thrust  augmentation 
of  the  XFV-12A  ejector  can  ba  increased  from  1.45  to  1.64  by  the  addition 
of  lobes  to  the  baseline  nozzle,  and  a corresponding  increase  of  throat 
width. 


INTRODUCTION 


The  static  thrust  of  turbojet  engines  can  be  significantly  increased 
by  diverting  the  exhaust  flow  through  an  ejector  pump.  According  to  the 
laws  of  momentum  and  energy  conservation,  greatest  thrust  is  obtained 
from  a given  energy  input  by  accelerating  a large  mass  of  air  to  a low 
exhaust  velocity.  Within  an  ejector,  thrust  is  increased  by  transferring 
the  kinetic  energy  of  the  engine  exhaust  stream  to  a larger  mass  of  air 
drawn  from  the  atmosphere.  The  ejector  duct  experiences  a reaction  force 
which  is  equal  but  opposite  to  the  momentum  change  of  the  accelerated 
stream.  Details  of  this  process  have  been  discussed  by  Bevilaqua.^ 

The  mechanism  of  this  energy  transfer  is  the  turbulent  mixing  of 
the  two  streams.  Thus,  increases  in  ejector  thrust  augmentation  can  be 
obtained  by  increasing  the  turbulent  mixing  rate.  Appreciable  increases 
in  mixing  and  augmentation  have  been  achieved  with  the  so-called  hyper- 
mixing2,3.  The  alternating  exit  of  the  hypermixing  nozzle  serves  to 
introduce  a row  of  streamwise  vortices  into  a plane  jet  (Figure  1). 


Figure  1.  Hypermixing  Nozzle  Exit 
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These  vortices  serve  to  accelerate  the  turbulent  mixing  and  thus 
to  entrain  additional  fluid  into  the  ejector.  Because  of  their  favorable 
entrainment  characteristics,  hypermixing  nozzles  in  the  centerbody  were 
used  in  the  XFV-12A  airplane.  Figure  2 illustrates  how  the  ejector 
components  fold  into  the  shape  of  a wing  for  forward  flight.  This  study 
is  concerned  with  the  development  of  the  centerbody  nozzle. 


CENTERBODY 

NOZZLE 


Figure  2.  Deflected  Ejector  wing 

The  specific  objective  of  this  investigation  was  to  develop  and 
demonstrate  centerbody  nozzles  that  would  provide  increased  augmentation, 
exceeding  the  peak  value  of  1.45  obtained  for  the  hypermix  configuration. 

In  the  following  sections,  a "master  plan"  is  presented  including 
a brief  description  of  the  3-dimensional,  turbulent  kinetic  energy  program 
(3D-TKE)  used  to  calculate  the  viscous  mixing.  Analytical  and  experimental 
results  for  both  symmetric  and  asymmetric,  centerbody  nozzle  configurations 
are  discussed. 


NOZZLE  DEVELOPMENT  PLAN 


The  overall  program  consisted  of  three  phases. 
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Phase  I To  Improve  Analytical  Capability  for  Calculating  3-Dimensional, 
Turbulent  Flow  by: 

•'  Identifying  inlet  turbulent  kinetic  properties  for 
augmenters  (k,  e) 

• Inlet  grid  generation  program 

• Simplification  for  incompressible  flow 

Phase  II  To  Compare  Analyses  with  Experiment  for  Hypermixing  Nozzles 
to  Update  Analytical  Technique 

Phase  III  To  Develop  and  Demonstrate  Centerbody  Nozzles  that  Increase 
Augmentation 

ANALYTICAL  MODEL  AND  NUMERICAL  APPROACH 


Governing  Equations 

The  program  equations,  in  cartesian  coordinates,  include: 
Equation  of  State 
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where  U,  V,  and  W are  the  time  averaged  components.  These  equations  are 
regarded  as  parabolic  in  the  longitudinal  coordinate,  X,  and  elliptic 
in  the  transverse  coordinates;  Y and  Z.  A more  complete  description  of 
the  program  is  given  by  DeJoode  and  Patankar^. 

The  turbulent  shear  stresses,  r^.j,  are  expressed  in  terms  of  a turbulent 
viscosity,  y^.,  and  velocity  gradients.  The  two-equation  turbulence  model  of 

'<  c 

Launder  and  Spalding  expresses  the  turbulent  viscosity  in  terms  of  two 
parameters,  K and  e,  for  which  two  differential  equations  are  solved.  The 
expression  for  is: 


(2) 

where  c^  is  an  empirical  constant  and  K and  e are  obtained  from  the 
solution  of: 


G is  the  rate  of  generation  of  K by  the  action  of  velocity  gradients. 
Since  the  only  significant  gradients  are  3U/3Y  and  3U/3Z,  the  expression 
for  G becomes 


(4) 
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Values  of  the 

constants  (Reference  5)  include: 

• c 

V 

"1 

^^2 

‘'K 

0.09 

1.44 

1.92 

1.0 

The  above  equations  were  put  in  finite  difference  form.  From  known 
conditions  at  an  upstream  cross  section,  X,  the  flow  field  at  the  down- 
stream cross  section,  X + aX  is  computed.  This  streamwise  marching 
process  is  continued  until  the  domain  of  interest  has  been  covered. 

Boundary  Conditions 

Figure  3 presents  a typical  ejector  configuration,  containing  a hyper- 
mix centerbody  nozzle.  Other  types  of  centerbody  nozzle  designs  can  also 
be  analyzed.  Needed  geometric  parameters  include  the  throat  to  'Inlet 
primary  flow  area,  A2/A0,  the  ejector  exit  to  throat  area,  A3/A2, 
the  diffuser  length,  L,  the  Coanda  surface  shape,  and  the  flow  split  or 
division  of  the  primary  flow  between  the  centerbody  nozzle  and  Coanda  jets. 


Figure  3.  Typical  Ejector  Wing  Configuration 


829 


The  computational  boundaries  are  outlined  in  Figure  4 for  the  hypermix 
type  configuration.  Symmetry  planes  are  used  as  computational  boundaries 
in  the  spanwise  direction  because  most  nozzle  designs  are  "periodic"  along 
the  span. 
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Figure  4.  Computational  Boundaries 
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Initial  Conditions 

Primary  jet  velocities  were  calculated  from  the  conventional  isen- 
tropic  relation  using  a velocity  coefficient  of  0.925  (reference  6).  The 
mean  static  pressure  along  with  the  ambient  pressures  are  used  in  Bernoulli's 
equation  to  calculate  the  inlet  secondary  stream  velocities.  The  inlet 
static  pressure  is  obtained  in  the  iterative  solution  of  the  viscous  flow 
• field  and  is  discussed  in  the  "Closure  Scheme".  Values  of  the  initial 
turbulence  conditions  (K,  e)  are  discussed  in  a later  section. 

Inlet  Grid  Generator 


At  the  inlet  computational  plane  of  the  ejector,  it  is  necessary 
that  the  flow  area  be  subdivided  into  a grid  of  approximately  1500  to 
3000  control  areas  (volumes)  with  corresponding  velocities  (primary  or 
secondary).  To  reduce  the  usual  20  hour  setup  time  to  an  hour  or  less 
and  also  eliminate  the  arbitrariness  of  the  subdivision,  a computer 
program  was  written  to  generate  the  grid  where  only  the  ejector  and 
noz2le  dimensions  are  the  primary  input.  A grid  for  the  hypermix  con- 
figuration is  shown  in  Figure  5 and  contains  approximately  3000  control 
volumes.  At  the  upper  and  lower  surfaces,  the  grid  appears  to  be 
shaded  because  in  the  Coanda  jet  region,  a much  finer  grid  is  used. 

Closure  Scheme 


Closure  is  obtained  by  iterating  on  the  inlet  pressure  until  the 
calculated  pressure  at  the  end  of  the  exhaust  plume  is  ambient.  The 
plume  exit  pressure  is  dependent  upon  the  curvature  of  the  jet  sheet 
leaving  the  trailing  edge  of  the  ejector  shroud  as  well  as  the  plume 
length,  both  determined  by  the  difference  between  the  ambient  pressure 
and  the  ejector  exit  pressure.  This  "jet  flap"  effect  is  discussed 
in  more  detail  in  reference  4. 

The  thrust  augmentation  ratio,  <}.,  is  defined  to  be  the  ratio  of 
the  ejector  stream  thrust  to  the  isentropic  thrust  obtained  by  expanding 
the  same  mass  of  primary  fluid  to  atmospheric  pressure.  The  thrust  of 
the  ejector  is  evaluated  by  integrating  the  exit  momentum  flux  and 
pressure  force. 


pU^  dY  dZ  - (P^  - P3)A3 


in  which  P^  and  A,  are  the  static  pressure  and  area  at  the  exit.  (()  is 
defined  as: 
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Simplification  for  Incompressible  Flow 

To  reduce  computer  time  the  analysis  was  conducted  using  the  incom- 
pressible option  in  the  program.  With  this,  density  is  assumed  constant 
and  the  state  and  energy  equations  are  bypassed.  To  verify  that  this  was 
a valid  approach,  both  compressible  and  incompressible  computer  runs  were 
made  for  a hypermix  configuration  at  a pressure  ratio  of  2.2  and  primary 
gas  temperatures  of  80® F (usual  test  conditions).  There  was  little 
difference  in  the  ejector  calculated  exit  veiiK-'tj  profiles.  This 
afforded  a saving  of  about  30. percent  or  10  m'lnutes  ot  CPU  time  on  the 
CDC-176  computer. 


RESULTS  AND  DISCUSSION 


The  following  sections  describe  how  the  hypermix  test  results  inter- 
acted with  the  computer  program  development  (Phase  II)  so  that  it  could  be 
used  to  identify  higher  performance  symmetric  and/or  asymmetric  nozzle 
configurations  (Phase  III). 

Hypermixing  Nozzles 

To  improve  the  analytical  design  methods,  angularity  measurements 
were  made  of  the  secondary  flow  at  the  ejector  inlet.  These  data  were 
generalized  and  incorporated  into  the  program.  Initial  values  of  the 
kinetic  energy  for  use  at  the  inlet  computational  plane  were  obtained 
from  hot  film  measurements. 

Figure  6 presents  the  calculated,  2-dimensional  mainstream  velocity 
distribution  along  with  the  calculated  and  measured  profiles  for  two 
lateral  positions,  corresponding  to  the  middle  of  the  hypermix  element 
and  between  adjacent  elements.  It  should  be  noted  that  the  Coanda  velocities 
are  not  in  good  agreemenx:  since  Coanda  curvature  effects  have  not  been 
included  in  the  computer  program.  In  the  2-dimensional  velocity  distribution, 
the  nearby  Coanda  velocities  were  purposely  omitted  from  the  plot  in  order 
that  the  intermediate  velocity  distribution  could  be  shown. 

Calculated  <}i's  are  compared  with  measured  values  in  Figure  7.  Analytical 
results  were  normalized  at  a of  1.45  at  7°  since  the  computer  program  pre- 
dicts a A({i.  The  7®  reference  point  corresponds  to  the  <{.  measured  for  the 
7®  hypermixing  nozzle  configuration,  used  in  the  XFV-12A  airplane.  Results 
are  for  an  A3M2  = A2/A0  = 16.2,  and  an  L/D  = 1.8,  where  D is  the 

throat  width.  Up  to  about  20®,  the  predicted  trend  agrees  with  measured 
values.  Beyond  20°,  calculated  ({.'s  decreased  while  the  measured  values 
remained  essentially  constant. 
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SECTION 


Figure  6.  Comparison  of  Calculated  and  Measured  Velocity  Profiles  at  Ejector  Exit  for  22.5°  Hypermix  Nozzles 
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Figure  7.  Comparison  of  Measured  and  Computed  Augmentation 

• With  the  increasing  hypermixing  angle,  the  swirling  or  vortex  action 
should  increase  and  therefore  additional  swirl  terms  (originally  deleted 
to  reduce  computer  time)  should  be  included  in  the  turbulence  generation 
term,  G (Equation  4).  With  the  addition  of  all  swirl  terms  in  G for 
the  28°  case,  the  calculated  increased  by  only  0.01.  The  lack  of 
agreement  beyond  about  20°  is  attributed  to  a limitation  in  the  manner 
used  to  solve  the  flow  equations,  in  that  the  V and  W velocity  components 
should  be  much  smaller  than  the  mainstream  velocity,  U.  At  high  swirl 
angles  this  may  not  be  true. 

Since  good  agreement  of  calculated  and  measured  velocities  ere 
obtained  along  with  the  general  ^ trends,  the  computer  program  was  used 
in  the  design  and  analyses  of  other  nozzle  configurations,  discussed  in  ■ 
the  following  sections. 


Asymmetric  Nozzles 


This  type  nozzle  was  developed  in  an  effort  to  capitalize  on  the 
generally  superior  entrainment  characteristics  of  the  symmetric  cross- 
slot nozzle  configuration  while  niainUining  packaging  limits  imposed  by 
supersonic  airfoil  contours.  The  design  combines  a series  of  aft  facing, 
spanwise  slot  nozzles  with  an  alternating  series  of  cross-slots  on  the 
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forward  side.  An  upstream  view  of  the  14  element,  asymmetric  centerbody 

nozzle  along  with  a side  v<ie.'  is  shown  in  Figure  8.  The  aspect  ratio  of  'O'. 

the  span  slot  is  5.5.  i 


Figure  8.  14  Element  Asymmetric  Nozzle 

The  computer  program  was  used  to  indicate  expected  4<j)trepds.  Tests 
were  then  used  to  confirm  the  predicted  trends.  For  example,  to  increase 
^ it  was  analytically  shown  that  for  this  type  nozzle,  it  is  necessary  to 
increase  the  throat  width,  A2,  to  prevent  the  merging  of  the  cross-slot 
and  Coanda  jets.  In  Figure  9 the  measured  maximum  <j)  for  the  5.4  aspect 
ratio  configuration,  wide  throat  was  1.53,  an  increase  of  0.02  over 
the  hypermixing  peak  value  (Figure  7). 
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Figure  9.  Effect  of  Span  Slot  Aspect  Ratio,  Asymmetric  Nozzles 

Further  analyses  showed  that  when  the  aspect  ratio  of  the  span  slot 
is  reduced  from  5.5  to  2.8,  the  (|)  should  increase  by  0.03.  Measurements 
(Figure  9)  showed  a gain  of  0.02.  The  peak  (j.  for  the  asymmetric  con- 
figuration was  1.55,  a a<|)  of  0.04  above  the  peak  for  the  hypermixing  nozzles. 

Calculated  and  measured  velocities  for  the  smaller  aspect  ratio  con- 
figuration are  shown  in  Figure  10.  Velocity  profiles  are  similar. 

Symmetric  Nozzles 

Based  on  the  previous  analytical  and  experimental  results  neither 
the  hypermixing  or  asymmetric  nozzles  provided  a (}>  of  the  order  of  1.65, 
the  target  value.  The  computer  program,  however,  has  been  p»"C 'en  useful 
in  their  design  and  evaluation.  To  obtain  reliable  results,  t'le  program 
does  require  an  accurate  description  of  the  flow  conditions  at  the  inlet 
computational  plane,  particularly  that  of  the  primary  flow.  Secondary 
flow  angularity  was  discussed  earlier.  Therefore,  the  study  of  the 
symmetric  nozzles  was  concerned  with  primary  flow  angularity  measurements, 
the  design  and  testing  of  a nozzle  to  capitalize  on  these  measurements, 
and  finally  the  effect  of  geometric  variations  on  (j>. 
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The  symmetric  nozzle  configurations  consist  of  either  a series  of 
symmetric  cross  slots  or  a series  of  cross  slots  with  intennediate 
span  slots.  Figure  11  shows  schematically  the  downstream  vortex  action 
for  symmetric  cross  slots,  in  the  absence  of  span  slots. 


Figure  11.  Symmetric  Cross  Slot  Flow  Pattern 

Figure  12  presents  a typical  cross  slot-span  slot  configuration. 

In  this  configuration  the  thickness  of  the  cross  slot  gap  is  shown  to 
vary  linearly  from  the  centerline.  In  other  configurations  this  thick- 
ness is  held  constant.  Geometric  design  parameters  are  shown  in  Figure  13. 


Cross  Slot  Primary  Flow  Discharge  Angle:  Angularity  measurements  were 
made.  Figure  14  presents  the  data  for  four  nozzle  configurations  having 
the  following  wedge  and  launch  angles:  (a)  0°,  21®;  (b)  21®,  21®;  (c)  11®, 

17.2®;  and  (d)  21®,  9®.  Schematics  of  the  corresponding  nozzles  are 
also  shown.  The  abscissa  is  the  nondimensional  distance  from  the  center- 
line  to  the  edge  of  the  nozzle.  Measurements  were  made  under  four  nozzles 
in  each  configuration.  There  is  some  scatter  due  to  slight  manufacturing 
differences  for  the  small  scale  model  augmenters,  evident  in  Figure  14(c). 

For  analytical  purposes,  the  profile  shapes  were  approximated  by  straight 
lines.  Calculated  A<j)'s  are  shown  in  Figure  15  for  a linearly  varying  and 
for  a constant  angular  distribution.  The  abscissa  is  the  average  of 
the  angle  at  S = 0 and  S = 1.  Analytically,  a nozzle  having  the  linearly 
varying  flow  angle  provides  a greater  Ai>.  To  verify  the  trend,  a cross 
slot  nozzle  v;as  designed  to  provide  this  type  profile  (Figure  16).  Figure 
17  presents  a comparison  of  the  measured  flow  angularity  with  the  desired 
linear  distribution.  The  results  demonstrate  that  a nozzle  can  be  designed 
to  provide  a desired  distribution.  For  this  configuration  the  predicted  A(|) 
was  0.30.  The  measured  value  was  0.24  (Figure  15). 

Bowtie  Ratio:  A "bowtie  nozzle"  shape  (cross  slot  nozzle  gap  varying, 
as  shown  in  Figure  12)  should  increase  entrainment  and,  hence,  (j>  by  placing 
more  primary  flow  into  the  "dog  bone"  vortices  at  the  cross  slot  tips 
(Figure  11).  Analytical  results  are  compared  with  experimental  results 
in  Figure  18.  Up  to  a bowtie  ratio  of  about  2,  (j.  does  increase  as  predicted. 

Span  Slot-Cross  Slot  Flow  Split:  Both  analytical  and  experimental  tests 
were  conducted  with  the  cross  slot-span  slot  configurations  to  determine 
the  effect  of  flow  split.  Flow  split  is  defined  as  the  ratio 
of  the  flow  through  the  span  slots  to  the  total  primary  flow.  Results 
are  compared  with  data  in  Figure  19.  Best  performance  is  obtained  when 
the  flow  to  the  span  slot  is  about  40  percent  or  less.  These  results 
are  understandable  since  the  downstream  jet  expansion  from  the  span  slot 
can  generate  destructive  interference  with  the  flow  from  the  cross  slot. 


Cross  Slot  Aspect  Ratio:  Since  the  results  in  the  previous  section 
showed  that  the  percent  primary  flow  to  the  cross  slot  nozzle  should  be 
increased,  the  next  logical  step  was  to  investigate  the  effect  of  cross 
slot  aspect  ratio.  Earlier  work  by  Peschke^  showed  that  the  entrainment 
for  2-dimensional  free  jets  is  increased  with  aspect  ratio.  Predicted 
results  for  augmenters  (Figure  20)  show  that  aij)  is  increased  with  cross 
slot  aspect  ratio;  however,  physical  width  constraints  must  be  considered. 


A(j> 


Figure  20.  Effect  of  Cross  Slot  Aspect  Ratio 

Symmetric  Nozzle  Overall  Performance:  By  combining  the  computer  program 
with  tests  of  the  inlet  flow  angularity,  geometric  variations,  etc.,  a 
'measured  peak  (|)  of  1.64  was  demonstrated  for  a cross  slot-span  slot 
configuration  (Figure  16). 


CONCLUSIONS 


The  following  results  were  obtained  from  a combined  analytical 
and  experimental  investigation  of  centerbody  nozzle  configurations. 

. By  combining  the  data  of  hypermixing  nozzle  configurations 
with  a 3-dimensional,  turbulent  kinetic  energy  computer  program, 
an  analytical  procedure  was  developed  for  predicting  the  flow 
field  in  ejectors  and  their  augmentation,  ({>. 

. The  measured  <i)  for  the  7°  hypermixing  nozzle,  used  in  the  XFV-12A 
airplane,  was  1.45.  It  was  analytically  and  experimentally  shown 
that  a peak  <}>  of  1.51  could  be  obtained  by  increasing  the  hyper- 
rnixiiig  angle  to  22". 


. Data  for  the  asymmetric  nozzle  showed  a peak  ^ of  1.55. 


. Although  larger  <|)'s  were  obtained  by  extending  the  width  of  the 
cross  slot,  optimization  studies  and  tests  showed  that  to  increase 
(|)  it  Wft>  necessary  to  increase  the  bowtie  ratio,  relative  flow  to 
the  cross  slot,  cross  slot  aspect  ratio,  and  to  provide  a cross 
slot,  exit  velocity  angularity  which  varied  linearly  with 
distance. 

. By  combining  the  computer  program  with  test  data  of  inlet 
flow  angularity,  geometric  variations,  etc.,  a measured  peak  (|> 
of  1.64  was  demonstrated  for  the  cross  slot-span  slot 
configuration  in  which  the  exit  angularity  of  the  cross  slot 
varied  linearly  with  distance. 
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